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y < M321/H 
Section A 


Answer no more than THREE questions from this section. 
All questions in this section have equal marks. 


Question 1 


Write down the function which satisfies Laplace's equation in a square 
and which takes the value 1 on all four sides. 


Show, by using arguments of symmetry, that if u(x, y) satisfies Laplace’s 
equation in a square, whose centre is at the origin, and takes the value 0 
on three sides of the square and 1 on the fourth, then u(0, 0) = 4. Explain 
why your approach is valid. : 


Hence, by finding the solution u(x, Y) throughout the square, verify that 


|a 


l 


o 
(= 1) 1 


Question 2 

Show that for a scalar function u with domain R? 
div(u grad u) = (grad u)? + uV2u, 

and using the Divergence Theorem, or otherwise, deduce that if 
(a) V?u = — 1 in a volume D, 


(b) at the bounding surface, S, of D either u = 0 on all’S or 
Gu/dn = O on all S, . 


fff, wr>o. 


Question 3 
Find the solution of the problem 


then 


2, 
yw) = 3 wax O<x<1, 


w(0) = o(1) =0, 
in terms of the eigenfunctions of L or otherwise. 
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Section A — continued 


Question 4 
The following problem, 


ht a? 
wtp’ O<x<10<y< 4, 


u(x, 0) = x O<x<il, 
u(0, y) = 0 O<y<}, 
ul, y) = 1 O<y<}t, 
F(x, 4) =0 O0<x<l, 


is to be solved by the finite-difference method using a uniform square 
mesh of side }. 


What are the numerical schemes when the derivative u/dy is replaced by 
(a) a forward-difference formula? 
(b) a central-difference formula? 

If the numerical scheme in case (b) is put in the form 
Au = b, 


where A is a matrix whose diagonal elements are 1 and b is a vector 
whose elements are known, can the simple iterative scheme given by 


uw”) = (Aue) +b 


be used to solve the equations? 


a“ 
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Section B 


Answer no more than TWO questions from this section. 
All questions in this section have equal marks. 


Question 5 

A long elastic string is moving with velocity v in the x-direction; it passes 
through two small tings fixed to the x-axis at the origin and x = I. By 
considering the equation of motion in a frame of reference moving with 


the string, show that the transverse motion of the string, for 0 < x <l, 
is described by the equations 


( =) ay 2v ay 


1 @y 
l-a 3 


ax? ~ 8 dea oF EO» 
WO, 1) = (1,2) = 0, 


where c is the speed at which small disturbances propagate along the 
stationary string. 


If v < c, show that the natural frequencies of vibration are 


es nnc(l — v?/c?) 


and that at each of these frequencies the corresponding solution will be 
Proportional to the real part of 


nD = epi (Fife + Diete- 98) 


- exp (F ore - 1) {x -(c + na). 


Interpret this equation in terms of wave motion. > 


(No marks will be given for simply substituting y, into the differential 
equation.) a 


Qusstfon 6 
For 0 < x < x show that if a is not an integer then 


kad 
_ sinan | sinar af 12 1 
wawe he > (= 0) (it Ay) com, 
a 


and deduce that 


1 1 1 
rom- ie > (+5): (1) 


aad 


By showing that 
lim (cot me — 2) de = in 


a0 e 


sin xz 
RZ 


» 


deduce that for |z| <1, 
sin xz 


a ea 9(1-2)(1-2)...(r-2). 


(You may assume that the series in equation (1) is uniformly convergent 
on the interval [— p, p) for any pe (0, 1).) 
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Section B — continued 


Question 7 
Write the equation 


2, 
Sa 2H 4 u= Fo) 


in self-adjoint form. 
Show that the influence function is given by 
R(x, €) = e“ GF sin(x — €), 
and obtain the solution satisfying u(0) = a, u'(0) = b in terms of the in- 
fluence function. 


Deduce that if F is bounded on (0, œ) then all solutions Satisfying the 
given initial conditions are bounded on [0, 0). Suggest an argument to 
show that if in addition F(x) ~ 0 as x ~ œ then u(x) ~ 0 in this limit. 


Question 8 
Classify the equation 
3u eu Ou 
-lP — x?" a a-1 T5 + 

(n Dirr x" a nx? Fe œ@NeERx R*, nél. 
Reduce it to standard (canonical) form and write down the general 
solution. Hence show that the solution which satisfies the initial con- 
ditions 


MOR), Eea) xer, ; 


ux, 1) = HRW.) + Fwy) ES Dae 
where B = 1 — n, wg = (x° + QP. 
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Section B — continued 


Question 9 
For the diffusion equation 


au aU 
- Ot ax? 


Lees’ scheme is defined as 
(1 +r) Up, get — Mpa, 4+1 = (1 — F) uy + Tups, q> 


wherer = k/h? isconstantand k and h are the mesh spacings in the /- and x- 
directions respectively. Using von Neumann's method, find the conditions 
under which Lees’ scheme is stable. 


(N.B. If z,, zz and € are complex numbers such that € = z,/z, then 
Izu] = 11 |zal.) 


What conditions must be imposed on / and k in order that the resulting 
numerical solution is convergent to the solution of the differential equa- 
tion? 


mY 


